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Abstract 

The vacuum energy of two CP{1) sohtons on a torus is computed numerically. 
A numerical technique for the zeta-function regularisation is proposed to re- 
move the divergence of the vacuum energy. After performing the numerical 
regularisation, we observe the effect of the vacuum energy on the two-soliton 
configuration. 
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I. INTRODUCTION 



Topological solitons arise as solutions of classical field equations and corresdond to points 
where the classical energy attains degenerate local minima. The classical energy of a topo- 
logical soliton typically depends on an integer which can be interpreted as the number of 
solitons, and in some cases it is independent of the parameters which correspond to the 
positions and the sizes of individual solitons. Therefore the quantum corrections can play a 
dominant role in the interactions between solitons. 

Many techniques have been used for the numerical computation of the vacuum energies of 
single solitons . These methods rely heavily on the spherical symmetry which is usually 
implicit in the single soliton solutions. In this paper, we propose a new technique which 
does not rely on any symmetry and can therefore be applied to multiple soliton solutions. 

We will evaluate the vacuum energy for the CP{1) model in d = 1 + 2 dimensions. 
Previous work on single solitons in the CP{1) model has used an approach based on heat 
kernel coefficients and phase shifts [0. It was shown that the solitons are unstable to 
collapse due to the quantum corrections. The CP(1) sigma model also exhibits multi- 
soliton solutions. We focus on the charge-2 case and study their interaction due to quantum 
corrections. 

The quantisation of CP(1) solitons can be performed by following the standard pertur- 
bation techniques invented by Schwinger and developed in the 1970s by several authors 
(for a review, see 0). The regularisation will be performed by the zeta- function technique, 
which we convert into a form that can be evaluated numerically. We examine how the finite 
one-loop energy depends on parameters such as the separation of the two solitons and their 
width. 



The CP{1) sigma model consists of a single complex scalar field taking a value in the 
one dimensional complex projective plane CP{1). The action of the CP{1) sigma model is 
given by 



where fi is the spacetime index running over 0, 1,2 and |mP = uu. The complex projective 
plane being mapped to a sphere stereographically, the CP{1) model is equivalent to the 
0(3) model. Let u and {4>°'}, a = 1, 2, 3 be coordinates of the CP{1) and 0(3) respectively. 
The stereographic coordinates are given by 



II. CP{1) SIGMA MODEL IN (2 + 1) DIMENSIONS 




(1) 




(2) 



with u = Ui + iu2 and (0")^ = 1. Expressing the action (|l|) in terms of (f)"", one can recover 
the 0(3) sigma model action 




(3) 
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Topological soliton solutions for the 0(3), or equivalently CP(1) model, were discovered 
by Belavin et al. §]. To be precise, their solutions are instantons in two-dimensional Eu- 
clidean spacetime, but at a classical level, they are the same as solitons in three-dimensional 
Minkowski spacetime. 

The static energy of the 0(3) model is given by 

E=\jd?x {d.rf (4) 
where i takes the values 1, 2. The finite energy condition requires a boundary condition 

dicf)" ^0 as |x| ^ oo . (5) 
Without loss of generality, we can define the asymptotic value of 0" as 

(0")2 ^1 as |x| ^ oo (6) 

which satisfies (|). This boundary condition compactifies the space into a sphere. Since the 
field space is also a sphere, the homotopy group of the fields is vr2(S'^) = Z and hence there 
are soliton solutions. The topological charge is given by 

Q = ^ j £x e.,ea6c0"(9./)(9,0'=) . (7) 



57r 

From the obvious identity 

{d.r - e^,eafec0'9,-0^)2 > , (8) 

one can derive 

{d,rf>e,,,abcr{d,(t)'){dj(t)'). (9) 

This implies 

E > 2tt\Q\ . (10) 

The soliton solutions attain minimum energy in each topological sector and hence saturate 
the equality. Thus they satisfy the first order differential equation 

- eijeabcct^'djct)' = . (11) 

This equation turns out to be equivalent to the Cauchy-Riemann conditions in the OP(l) 
version 

dui _ du2 du2 _ _dui^ 
dx dy ' dx dy 

The general solution is given by the analytic function 

mz) = 7 ^ 7 r 13) 
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where z = x + iy^ and Oj, {i = 1, ■ ■ - n) and bj, {j = 1, ■ ■ - m) are the complex parameters 
characterising the position and size of sohtons respectively. Note that the only singularities 
the field u can have are isolated poles, since they merely correspond to the north pole of a 
sphere in the target space of the 0(3) model. 

The degree of the mapping is equal to the topological charge of the solution (p!3|), which 
is also equal to the number of solutions when expressing z in terms of mq, 

Q = max(m, . (14) 

u corresponds to the anti-soliton solution and gives an opposite charge of u. 



III. CHARGE-2 CP(1) SOLITONS ON THE TORUS 

For the sake of numerical work, we consider the CP(1) model on the torus [Q. This 
requires a slight modification of the previous soliton solutions. We shall assume that the 
field u{z) on the torus satisfies the periodic boundary condition 

u{z + 1 + i) = u{z) (15) 

The fact that the only singularities of Uq{z) are poles means Uo{z) is a meromorphic function. 
Then Uo{z) can be represented as an elliptic function. 

In the charge-2 case, the solution is represented by the Weierstrass function (see appendix 
For simplicity, we shall assume that the two solitons are symmetric in their size and 
location. Then there are only two complex parameters required which are the degrees of 
freedom of the size and separation. Thus we may write 

1 

with complex parameters a and p. 

Let the separation of the solitons (located at the poles) and the width be 2e and w 
respectively and restrict to real parameters. For (0), one obtains 

2 

P = -P(^e), a = — - . (17) 

wp'[te) 

In terms of e and w, the solution (ITBD becomes 



, , wp he) 

uJz) = , , : \ . 18) 

The energy density for this solution with e = 0.3 and w = 0.5 is plotted in fig.|l|. 
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FIGURES 




0.4 



FIG. 1. Energy density of charge 2-soliton for e = 0.3 and w = 0.5. 



IV. VACUUM ENERGY OF CHARGE-2 CP{1) SOLITONS 
The static energy functional for the CP(1) model on the torus is given by 



E\u] 



d X 



\d.iu\ 


2 


(1 + 


u 


2^2 



The analytic expansion at Uq is 



E[u] = Eo+ / (fxG^^f^ + - 



(19) 



(20) 



where ^(z, z) = u{z, z) — Uq{z). The 'metric' G and the fluctuation operator Aj are defined 
by 



G 



1 



and 



_y2 ^ 4wo(g^Mo) ^^' 

1 + |moP 



(21) 



(22) 



respectively. The computation of the vacuum energy depends on solving the eigenvalue 
problem 



(23) 



Zero eigenvalues correspond to changes in the parameters, or moduli, of the soliton solu- 
tion. The quantisation of these parameters is treated separately. Expanding ^ in terms of 
the remaining eigenfunctions gives an infinite sequence of (degeneracy two) oscillators with 
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vacuum energy A^^. The total vacuum energy can be regularised using the zeta- function 
scheme |TB|,ITB|, 



E = Eo + ({-l) (24) 
where the sum over non-zero eigenvalues 

oo 

C(^) = Y.K'- (25) 

n=l 

gives the generalised Riemann zeta function for s > 1. The value at s = —1/2 is uniquely 
determined by analytic continuation. 

V. NUMERICAL REGULARISATION 

We will now show how the zeta function regularisation scheme, which depends on analytic 
continuation, can be converted into a numerical subtraction scheme. Suppose that on average 
the eigenvalues approach the form 

Xn = an + b . (26) 
Using (j2^), we can define the corresponding heat kernel 

oo 

K{t) = J2 e"^"* • (27) 

n=l 

The behiour of the heat kernel as t ^ gives information about the eigenvalues at large 
values of n. We have 

where bn{x) are the Bernoulli polynomials and the first three terms are given by 

1 1 

bo{x) = l , bi{x) = X - - , b2{x) = x^ - X + - . (29) 

For small t, the leading terms are 

Comparing the order of t in ( p^ and (pO|), one can deduce 



where Bq and Bi are the heat kernel coefficients for the operator ([l8| ) calculated in appendix 
B, 5o = 1/ (47r) and Bi = 2. Thus, a = Air and b = — Gvr, with the result that 
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A„ = 4?™ - 67r . (32) 
We assign a zeta function to this eigenvalue 



C(^) = J:K'=(^) Us, -) (33) 



n=0 

where Ch{s,\) is the generahsed zeta function defined by 



CH(s,a) = ^(r^ + a)-^ (34) 

n=0 



Hurwitz has given a proof of the following formula 

2V{l-s 



CH{s,a) 



(27r)i-^ 
"2 J 



STi 



sin ( — j ^ n'* ^ cos(2 
^ ^ n=l 



+ cos I — ) ^ n*"""^ sin(2n'7ra) 

n=l 



Using this expression, one obtains 



, /l\"2r(l-s)~ ,1, . /S7r\ 



Setting s = —1/2 gives 



1 



oo (' 1 \n+l 



We can obtain a finite sum by subtracting the divergent terms, 

N 

N^Qo ■ 



iV 

Geg(s) = jim E(V-V) (37) 

n=l 

where all the zero modes are removed from the eigenvalues. Then one can write the true 
zeta function as 

C{s) = Crcs{s)+as). (38) 

Analytic continuation now implies 

The above expression allows us to evaluate the one-loop energy of the solitons numerically, 
assuming that (37) can be evaluated numerically. In practice, there are problems with the 
limit ^ oo at s = — 1/2 because the residuals A„ — A„ do not vanish as ^ oo. However, 
we shall see in later sections how the difference does vanish 'on average', allowing a result 
to be obtained. 
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VI. NUMERICAL COMPUTATION 



We use the Jacobi theta functions as a representation of the Weierstrass p function (see 
appendix ^ 



p{z) = n 



Thus the charge 2-sohton solution (|18D can be written as 



(40) 



Uo{z) = ZTTW — ^ — - . (41) 

In the following numerical computation, we shall use (^TJ) as a soliton solution. 

We shall evaluate the approximate eigenvalues uj in (^31) by the Rayleigh-Ritz variational 
method. This method introduces two principles. Firstly, the eigenfunctions are stationary 
configurations of the functional 

_ It, G(A,i 



Thus the eigenfunctions satisfy 

5^ 



. (43) 



And secondly, trial functions which give stationary values of E provide upper bounds for 
the eigenvalues (Hylleraas-Undheim theorem). 

The trial functions can be constructed as follows. We take N linearly independent 
functions 0„ parametrised by N variational parameters c„ and construct the trial function 
as their superposition 

N 

e = E Cn0n . (44) 
n=l 

Inserting into (H2) one obtains the functional E a.s a. function of the variational parameters 



Sm,?i=l CmCnBn 



E[ci,C2,---cn\ = ^ ^ ^ „ (45) 



^mn = / (fxG(t)m^f(t)n (46) 



where we defined 



Bmn = d X G(j)m4>n ■ (47) 

Then, from the minimum principle, the upper bound of the nth eigenvalue is given by the 
stationary point of E[ci, C2, • ■ ■ cn] with respect to c„, i.e. E satisfies 
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dE[ci,C2, ■■■cn] 



. 



(48) 



Thus we obtain linear homogeneous equations for each upper bound of (n = 1, ■ ■ ■ A^) 



N 







(49) 



m=l 



and the problem is reduced to the nth-degree secular equation of the N x N matrix 

det(A-EB) = 0. (50) 

Increasing in the number of basis states gives a lower upper bound of the exact eigenvalue 
in each mode. 

We take the basis of trial functions 0^ to be plane waves 



exp |^(A;z + kz)^ 



(51) 



where k = j{n+im) = 2Tx{n+im) and k = (n,m), n,m = —N, — ■ ■ ■ , 0, ■ ■ ■ , A^ — 1, A^. 
n and m are the mode numbers in the x and y directions respectively. With this basis, the 
functional fH5|) becomes 



Bk'k 



where 



k'k 



Bk'k 



dzdz- 



1 



2^2 



k'k exp 



-{{k' - k)z + {k' - k)z] 



dzdz- 



2^2 



exp 



■{{k' - k)z + {k' - k)z} 



(52) 

(53) 
(54) 



Therefore the upper bound of the spectrum can be computed from the secular equation of 
an (2A^ + 1) x (2A^ + 1) matrix 



det(Ap^ - EB 



k'k) 







(55) 



The spectrum will be more accurate as A^ increases and in the limit A^ — * cxd, it is exact. 

To solve (|55|) numerically, we simply used the LAPACK (Linear Algebra Package) which 
provides routines for solving systems of simultaneous linear equations, least-square solutions 
of linear systems of equations, eigenvalue problems, and singular value problems. The file 
actually used is chegv.f for computing all eigenvalues of an Hermitian-definite generalised 
eigenproblem. 
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VII. NUMERICAL RESULTS 

We have plotted the first 200 eigenvalues in fig.(|^) as a function of the mode number 
n = [k, k) for a typical background solution. The agreement with the asymptotic formula 
(^) is quite striking. Fig.@ is a plot of the regularised zeta function as a function of the 
number of modes up to which we took the sum of the eigenvalues. The fluctuations in 
the data reflect the fact that the eigenvalues deviate from the asymtotic formula. Although 
we can easily improve the accuracy of the individual eigenvalues, these fluctuations have a 
significant influence in preventing the convergence of the zeta-function in all of the cases 
which we have examined. This problem has been seen before in numerical calculations of 
zeta-functions [|l^. Averaging over has been proposed as a solution to the problem and 
we use this in our calculation. 

Xn 
2500 T 
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0- 

50 100 150 200 

FIG. 2. nth-eigenvalue A„ as a function of mode number n. 
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FIG. 3. Regulated zeta function as a function of the maximum mode number N. 
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Fig.(^) is a contour plot of the average value of the zeta function on the parameter space 
of e and w. The data shows that the more spiky and closer two solitons are, the lower the 
vacuum energy is. Thus the classically stable two solitons become unstable to collapse and 
merger when the one-loop quantum correction is taken into consideration. 
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FIG. 4. One loop energy of two separated solitons as a function of e and w. 

As a special case of a charge-2 soliton, we also examined the one-loop energy for the 
solution 



w 



(56) 



This solution exhibits a ring-shape energy configuration as is shown in fig.(^. The numerical 
one-loop energy as a function of w is given in fig.(^. The plots denote the numerical values 
and the line is its interpolated function which is found to be 

0.269v^ . (57) 

The rotational symmetry implies that this is also a case that can be analysed using the 
phase shift technique 0. Preliminary results suggest an identical functional form for the 



vacuum energy, but the phase shifts give a different overall scale |]T7 . 
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FIG. 5. Energy density of a ring for w = 1.0. 
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FIG. 6. One loop energy of a ring as a function of w. 



VIII. CONCLUSION 



In this paper, we have presented a numerical method for the computation of finite one- 
loop energies for solitons and studied its effect on the interaction of two CP(1) sohtons. The 
method does not require any symmetry to be present in the problem, but it is expensive 
in computer time and relies on the averaging out of residual terms in the eigenvalue spec- 
trum. Further analysis of these eigenvalue residuals, which should be related to fundamental 
periods of the torus, would be useful and might improve the technique. 

The results are very satisfactory and show that initially separated solitons are unstable 
to collapse and merger. Furthermore, the results are in qualitative agreement with the fast 
and reliable phase shift method [0] for a ring configuration, which has circular symmetry. 

Interestingly, the vacuum energy is not the only effect on the dynamics of two interacting 
solitons. There are additional forces arising from the reduction of the quantum field theory 
to the CP{1) moduli space |T^. In this example, an additional term is given by the Ricci 
tensor of the moduli space, and gives a repulsive force between the two solitons which is 
larger than the effect of the vacuum energy at very small separations IT? . 
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In principle, our methods are applicable to other solitons such as skyrmions, although 
moving from two to three dimensions would hughely increase the computer time. The 
one- loop correction to the calculation of the skyrmion mass is significant. For physical 
values of the pion decay constant, the classical skyrmion has a mass about 50% larger 



than a nucleon mass. However, according to Moussallam and Kalafatis ||Tl|, including one 
loop energies within the framework of chiral perturbation theory, the nucleon mass can be 
correctly predicted to within 20%. Presumably, quantum effects within chiral perturbation 
theory would also have an effect on the force between two nucleons which could be calculated 
by the method presented here. 



APPENDIX A: WEIERSTRASS ELLIPTIC FUNCTIONS 



An elliptic function / is a function such that it is doubly periodic with two primitive 
periods 2uji^2uj2 whose ratio is not real, i.e. 



where m, n are integers and 



f{z + 2muJi + 2nuj2) = f{z) 



Imf^^^l^O. 



(Al) 



(A2) 



Thus for an elliptic function /(-z), the 2-plane can be partitioned into parallelograms with 
vertices Zq + 2ma;i + 2nuj2- The only singularities of f{z) allowed in a period parallelogram 
are poles. 

The weierstrass function p is an even elliptic function of order 2 with one double pole at 
z = 2mu!i + 2nuj2 and defined by 



p{z) = p{-z) 



1 



+ E 



1 



1 



{z — muJi — {muJi + nuj2y 



(A3) 



where m,n takes all integers except for m = n = 0. The series ( |A3|) converges everywhere 
except at the poles. 

The function p{z) satisfies the differential equation 



4^p{zf - 92p{z) -93 = 4:{p{z) - ei){p{z) - e2){p{z) - 63) 



where g2 and are invariants and determine the periods ui and 002 as 

1 1 



(72 = 60^ 



^73 = 140^ 



Ci, 62 and 63 are then given by 

ei = p{ui), 62 = p{uJ2), 63 = p{uJ2) ■ 

with 



{muJi + nuj2y 



(A4) 



(A5) 



(A6) 
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92 93 
ei + 62 + 63 = 0, 6162 + 6263 + 6361 = - — , 616263 = — . 

The Weierstrass function can be represented in terms of Jacobi theta functions as 



(A7) 



p(z) = tj + 



TT 



J = 1,2,3 



where 



TTZ 



(A8) 



(A9) 



and the theta functions are defined by 

00 

Qx{z, q) = ei{z) = 2q-^ ^ sin{(2n + 1)^} 

n=0 

00 

e^iz, q) = e2{z) = 2q^ ^"^"^'^ cos{(2n + l)z} 

n=0 

00 

^3(^, q) = esiz) = l + 2j2q^" cos(2nz) 

n=l 

00 

9iiz,q) = e^iz) = l + 2^(-l)V'cos(2n2) , 

n=l 

with q = e*'^2/'^i^ g_^^ and 63 can also be expressed by the theta functions 



TT 



TT 



TT 



Ver^=— ^4(0), v/e^^=— ^2(0), v/e7^=— ^3(0) 



2cji 



2cJi 



2uJi 



(AlO) 



We adopt the theta-function representation for the weierstrass function since the theta 
functions converge rapidly for n and the series are periodic. For definiteness, we take the 
Lemniscatic case 



uji = —iuJ2 = - and 62 = 



(All) 



Then from ( |A1CI| ), we can determine numerical values of 61 and 63 as 



61 = -63 = n'etio) = — r 



Choosing j = 2 in (|A8|) , one obtains 



TT V4 



^2(0)^4(0)^3^ 



6.875 



(A12) 



Oiinz) 



(A13) 
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APPENDIX B: HEAT KERNEL COEFFICIENTS 



The small time expansion of a heat-kernel provides us a great deal of information about 
operators and their eigenvalues. This expansion was substantially developed in the 1960s, 
and more of the mathematical details can be found in the original work and reviews for 
example [|12|,|14 . 



Consider the eigenvalue problem on a manifold M 

A0 = X(f) (Bl) 

The field can have both spacetime, spinor and internal group indices, while A is the elliptic 
operator. The (integrated) heat-kernel of the operator A is defined by 

n 

We consider operators of the form 

-D'^ + X (B3) 

with a gauge covariant derivative D = V + A acting on fields which are associated with a 
representation of some given gauge group. For these operators, Gilkey has shown that the 
function K{t) has an asymptotic expansion 



oo 



n=0 

in (i-dimensional space. In general, the Bn coefficients depend on the operator, the geometry 
of the manifold and the boundary conditions satisfied by the fields 0. 

The coefficients Bq and Bi for the CP{1) model can be obtained by applying the general 
expressions. First, we rewrite the fiuctuation operator obtained in in a covariant form 

_a.a„ + ii^i^._2(fl.fl, + fl,fl.)-^^i!f^ (B5) 

where is a covariant derivative defined by 

D^u = d^u - —— — . (B6) 

1 + For 

Then the heat kernel coefficients for this operator are 

Bq = ^ I d^xii{l) (B7) 
47r J 

Clearly, for a torus of unit area, Bq = l/(47r). For a solution of the type 

Mz) = , (B9) 

one can show directly that 

Bi = 2. (BIO) 

The same result can be seen to hold for the two-soliton solution on the torus, because the 
integral (p^ ) is the index of the mapping given by u. 
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